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We have systematically investigated the energy resolution of a magnetic micro-calorimeter (MMC) for atomic 
and molecular projectiles at impact energies ranging from i? Ri 13 to 150 keV. For atoms we obtained absolute 
energy resolutions down to AE r; 120 eV and relative energy resolutions down to AE/E Ri 10“^. We also 
studied in detail the MMC energy-response function to molecular projectiles of up to mass 56 u. We have 
demonstrated the capability of identifying neutral fragmentation products of these molecules by calorimetric 
mass spectrometry. We have modeled the MMC energy-response function for molecular projectiles and 
conclude that backscattering is the dominant source of the energy spread at the impact energies investigated. 
We have successfully demonstrated the use of a detector absorber coating to suppress such spreads. We briefly 
outline the use of MMC detectors in experiments on gas-phase collision reactions with neutral products. Our 
findings are of general interest for mass spectrometric techniques, particularly for those desiring to make 
neutral-particle mass measurements. 
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I. MASS SPECTROMETRY FOR MULTI-KEV 
NEUTRAL PARTICLES 

Mass spectrometry is one of the most important 
analytical techniques used in both fundamental sci¬ 
ences and industry applications. The diverse fields 
of use include protein and genome characterization 
in molecular biologjd^, determination of fundamental 
constants^, atmosphere composition characterization for 
Eartb^ and other solar system bodies^, isotope dating^, 
drug discovery^, etc. Most molecular mass-spectrometry 
methods employ electromagnetic fields to determine 
mass-to-charge ratios of ionized particles^. Two major 
limitations arise from this approach. First, the mass- 
to-charge ratios may be ambiguous, e.g., singly charged 
monomers are indistinguishable from doubly charged 
dimers. Second, neutral fragments cannot be analyzed 
directly and a “neutral loss” must be accepted. To cir¬ 
cumvent these limitations, the neutrals studied can be 
ionized and then analyzed in a standard way using elec¬ 
tromagnetic fields^. Although “gentle” ionization meth¬ 
ods exist nowadays (e.g., electrospray ionizationi^) , the 
ionization often results in further fragmentation and in 
multiple ionization. Each of these processes can lead to a 
significant increase of complexity in the measured mass- 
to-charge spectruroii. Unfolding such data into a mass 
spectrum of the original sample may be complicatedi^ 
and often only qualitative information can be retrievedi^. 

Calorimetry represents an alternative approach to 
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measure the masses of neutral particles. Calorimetric 
detection can be based on ionization, such as in the sur¬ 
face barrier or scintillation detectors applied mainly in 
nuclear and particle physics^^, or on thermal effects, as 
we discuss below. In all of these approaches, the calori¬ 
metric signal is proportional to the energy set free by 
the stopping of an incident particle, mostly given by its 
initial kinetic energy. Masses can be deduced from the 
magnitude of the kinetic energy if the incident velocity 
is known. 

In molecular physics, calorimetric detection has been 
applied in particular to gas-phase collision studies using 
fast particle beams. For those studies many of the reac¬ 
tion products are neutral and their detection and mass 
measurement is essential in order to determine reaction 
cross sections and fragmentation branching ratios. A 
range of studies have been devoted to collision-induced 
dissociation (CID) of molecular and cluster ions in a 
fast beam, interacting with a stationary gas target (see, 
e.g. jiSiie for recent work). Extensive work has also been 
directed to the study of dissociative recombination (DR) 
of molecular ions with electrons, which for singly charged 
ions produces exclusively neutral fragments^i. DR is 
closely related to the process of electron capture dissoci¬ 
ation (ECD), which is widely applied in the field of mass 
spectrometry for sequencing multiply-charged biomolec- 
ular cation a^^d^ . 

The DR reaction on singly charged ions was studied in 
particular with fast ion beams in a storage ring, utilizing 
a merged electron-ion beams configuration^^ where the 
electron-ion interaction energy could be tuned to small 
(sub-eV) values. Considering this reaction on a model 
molecule ABCD+, the neutral fragmentation channels 
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may include, e.g., 

ABCD++ e” ^ A + BCD (1) 

^ AB + CD (2) 

^A + B + C + D. (3) 

In a fast ion beam, the parent ion velocity (correspond¬ 
ing to multi-keV or MeV ion beam energies) is typically 
much larger than the relative velocities that the frag¬ 
ments gain in the dissociation process (from the typical 
kinetic energy release of < 10 eV). As a result, the DR 
products arising from the beam overlap region have ve¬ 
locities close to that of the parent ions. After separat¬ 
ing the ion beam by means of magnetic or electrostatic 
fields the neutral DR products are collected by a kinetic- 
energy-sensitive detector. From the known parent ion ve¬ 
locity, and from the measured fragment kinetic energies, 
the fragment masses can in general be easily assigned. 
This method has been successfully used in combination 
with magnetic heavy ion storage rings^i^— operating at 
~ MeV ion beam energies and utilizing surface-barrier 
detectors (SBDs)^. The approach, however, is restricted 
to high beam energies as SBD detectors have a thin, 
but significant, insensitive surface layer which cannot be 
penetrated by low-energy particles^^. Moreover, even at 
MeV beam energies the energy resolution of SBDs is not 
sufficient to identify heavy neutral molecular fragments 
with 1 u resolution^!. 

Stored ion beams of lower energy have been attracting 
significant attention for some time now!!. One emphasis 
of these studies has been on complex molecular ions. In 
order to optimize the storage times of such ions and min¬ 
imize the influence of thermal radiation, ion beam traps 
and storage rings involving cryogenic cooling of the en¬ 
closing experimental vacuum chamber have recently been 
commissioned or are currently under development. In 
these electrostatic storage devices!!"— the experimental 
vacuum chamber is cooled to ^ 10 K. The beam energy 
for stored singly charged ions typically ranges up to some 
E ^ 30 keV. Higher ion beam energies, up to 300 keV for 
singly charged ions, will be used in the cryogenic storage 
ring (CSR) ^^'^^ at the Max-Planck-Institute for Nuclear 
Physics in Heidelberg, Germany. This ring is also being 
equipped for electron-ion merged beams experiments!! 
to allow for DR experiments on complex molecular ions. 
It is this CSR development that mainly motivates the 
present search for alternative calorimetric detection tech¬ 
niques for molecular particles. 

A promising approach for solving the needs for low ki¬ 
netic energy measurements is the use of cryo-detectors!!, 
such as Superconducting Tunnel Junction detectors 
(STJ) or micro-calorimeters. In these detectors the in¬ 
cident particle excites various types of energy states in 
the absorbing bulk. Depending on the detector type, 
this excitation can be transformed into electrical signal. 
For example, in micro-calorimeters the particle kinetic 
energy is transformed into heat and a sensitive tempera¬ 
ture monitoring then provides the measure of the kinetic 
energy deposited. For cryo-detectors, operation at low 


temperatures (< 1 K) is needed to limit thermal noise 
and thus resolve the measured low energies. The ad¬ 
vantages of cryo-detectors were recognized as early as 
1935 !!. Their broader use was first motivated by astro- 
physical applications, e.g., in dark matter searches!!, and 
for high-precision X-ray spectroscopy!!. 

In the 1990s the advantages of cryo-detectors were also 
recognized for classical mass spectrometry!!i!!. Aside 
from the benefits explained above, an additional motiva¬ 
tion for using cryo-detectors was their ability to resolve 
very slow massive molecules where standard microchan- 
nel plate detectors fail. Compared to microchannel 
plates, calorimetric detectors also offer a single-particle 
detection efficiency that could, in principle, reach val¬ 
ues very close to unity. On the other hand, the aspect 
ratio of the channels (typically about 0.7) limits the effi¬ 
ciency of microchannel plates for detecting single-particle 
impacts to significantly below unity for many applica¬ 
tions. Various types of cryo-detectors have been imple¬ 
mented in previous mass spectrometric studies, mainly 
in time-of-flight mass spectrometry (TOF-MS)!!"— . Re¬ 
cently, and independently of our own attempts, a cryo¬ 
genic STJ detector has been employed in a tandem 
mass spectrometer— following an approach similar to the 
calorimetric mass determination discussed above for gas- 
phase molecular collision and ion storage ring studies. 

When used for X-ray spectroscopy!! cryo-detectors 
can achieve an excellent relative energy resolution of 
AE/E Ri 10“"^ and an absolute resolution down to 
AE r: 1 eV. For detecting massive particles, however, the 
resolution degrades by orders of magnitude. In heavy-ion 
experiments performed at MeV and GeV energy, resolu¬ 
tions of not better than ~ 10“^ have been reached!!. For 
a-spectrometry studies researchers measured AE/E r 
5 X 10“^ dX E = 5.5 MeV—""—. However, energy res¬ 
olutions of only ^ 10“^ were reported for detecting 
atoms and molecules at energies of < 10 keV— ""— i!!!!!!!!, 
i.e., in the regime most relevant for mass spectrometry. 
Only a few studies have attempted to understand this 
issue. Andersen!! theoretically predicted the cryogenic 
calorimeter energy resolution for heavy ions. In his cal¬ 
culations he included various projectile energy loss pro¬ 
cesses in the detector absorber, such as backscattering 
and Frenkel-pair formation, and estimated the amount 
of energy which is not converted into heat and thus can¬ 
not be detected. The statistical nature of the losses 
then defines the broadening in the acquired energy spec¬ 
tra. Quantitatively, Andersen’s predictions for various 
detector types match to within an order of magnitude 
the energy resolutions measured later at MeV to GeV 
energies!!. More recently Horansky et al. have numeri¬ 
cally modeled the micro-calorimeter energy spectra from 
5.5 MeV a-particles!! and have concluded that the domi¬ 
nant spectral broadening is due to lattice damages. Their 
model underestimates the energy spread by a factor of 
^ 2 compared to the measurements. We are not aware 
of any study attempting to theoretically describe the en¬ 
ergy resolution of a cryo-detectors for massive particles 
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at < 100 keV energies. Similarly, to the best of our 
knowledge there are no systematic experimental studies 
of the energy resolution for cryo-detectors as a function 
of projectile energy and type for sub-MeV energies. In 
the energy range from ^ 10 keV to 5.5 MeV, even demon¬ 
stration measurements are missing. 

We therefore believe that the first step for understand¬ 
ing and eventually improving the cryo-detector energy 
resolution for mass spectrometry is a systematic study 
of the detector energy-response function for atomic and 
simple molecular projectiles at sub-MeV energies. To 
this end we have experimentally studied the energy res¬ 
olution of a micro-calorimeter detector for atomic and 
molecular ions and neutrals at energies from ^ 13 to 
150 keV. The rest of the paper is organized as follows: 
In Section m we briefly introduce the specific class of 
micro-calorimeters employed in our experiments, namely 
the metallic magnetic calorimeter (MMC). In Section Hill 
our experimental setup is described. In Section m we 
present the acquired energy spectra and a corresponding 
model resolving the various energy loss processes. Based 
on this model, we also demonstrate the importance of 
detector absorber materials. We summarize and present 
our future plans in Section IVl 

II. METALLIC MAGNETIC CALORIMETER (MMC) 
DETECTORS 

A. Detection principle 

Micro-calorimeters are cryo-detectors typically oper¬ 
ating at < 100 mK— . In this temperature regime the 
kinetic energy of a detected particle heats up the detec¬ 
tor absorber such that the increase of the absorber tem¬ 
perature is proportional to the original projectile kinetic 
energy. The response is linear to a high degree. In or¬ 
der to monitor the absorber temperature, various sensor 
types can be used. MMCs are operated with paramag¬ 
netic sensors whose magnetization at low temperatures 
is a monotonic function of temperature. The magneti¬ 
zation change is measured with a superconducting coil 
inductively coupled to a superconducting quantum inter¬ 
ference device (SQUID). A weak thermal link between 
the absorber and a cold thermal bath ensures that the 
absorber cools back to the bath temperature after the en¬ 
ergy readout. A detailed description of the MMC opera¬ 
tion principle and micro-fabrication can be found inS“— . 

B. Energy resolution 

In general the energy resolution of an MMC detector 
setup is given by the thermal noise in the detector ab¬ 
sorber and sensor, by the electronic noise in the signal 
processing system, and by the processes transforming the 
kinetic energy of the detected particle into absorbed heat. 
The thermal and electronic noise levels do not depend on 


the nature and energy of the detected particles and thus 
define the intrinsic energy resolution of the detector. For 
MMCs the intrinsic energy spread is typically <100 eV 
(full width at half maximum, FWHM) and under very 
favourable conditions energy spreads of < 1.6 eV have 
been achieved^i^. 

The kinetic energy transformation into heat is very 
efficient for X-rays. There the energy absorption pro¬ 
cess is dominated by the photoelectric effect, followed by 
the photoelectron efficiently heating the electron gas and 
subsequent energy transfer to the phonons^. Only a neg¬ 
ligible part of the energy is lost, e.g., by kinetic energy 
of secondary electrons escaping the surface of the ab¬ 
sorber. In contrast to detecting energetic photons, stop¬ 
ping atomic and molecular projectiles in the absorber 
may lead to additional processes, such as backscatter- 
ing of the projectiles, sputtering of absorber material, 
and absorber lattice defect formation. All of these non- 
thermal processes remove part of the projectile kinetic 
energy before it can be converted into heat in the ab¬ 
sorber and measured by the temperature sensor. More¬ 
over, the statistical nature of these processes results in a 
spread of the detected energies. Thus in general the en¬ 
ergy resolution of an MMC detector (and cryo-detectors 
in general) depends on the absorber material as well as 
on the nature and energy of the impacting particles. 


III. EXPERIMENT FOR DETERMINING MMC ENERGY 
RESOLUTION 

In order to experimentally determine the energy res¬ 
olution of an MMC detector we have measured kinetic 
energies of single ions from a nearly monoener- 
getic, mass-to-charge-selected ion beam. The accumu¬ 
lated single-particle calorimetric signals, after suitable 
calibration to yield E^, were used to derive energy spec¬ 
tra P„i(A„i, Aib, M) which, in turn, were used to deter¬ 
mine the corresponding energy resolutions. The mea¬ 
sured spectra are labeled by the nominal ion beam energy 
Aib and by the ion mass M (since the systems studied 
are simple enough to uniquely specify the ion type from 
its mass). 


A. Ion beam production 

Ion beams were generated using the high-voltage ion 
beam platform at the Max Planck Institute for Nuclear 
Physics in Heidelberg, Germany^. The various singly 
charged molecular cations were produced in a gas dis¬ 
charge from a Penning ion source. While still on the plat¬ 
form, they were accelerated in a first stage by ^ 13 kV 
and then filtered by their mass-to-charge ratio in a 90° 
dipole magnet. After passing through a set of apertures 
the ions were further accelerated in a second stage by 
the potential applied to the platform relative to ground. 
By adjusting the platform voltage, the final ion energies 
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were varied in the range from iJib ~ 13 to 150 keV. Af¬ 
ter leaving the platform the ion beam was focused and 
directed by an electrostatic quadrupole and a set of mag¬ 
netic deflectors onto the detector setup ^ 10 m from the 
platform exit. 

The nominal beam energy Ajb was determined as a 
sum of the two involved acceleration voltages used. The 
estimated systematic uncertainties in the Ajb scaling and 
offset are ±0.5% and ±300 eV, respectively. Here and 
throughout all uncertainties are quoted at an estimated 
la statistical confidence level. 

Some of the tests turned out to be sensitive to insta¬ 
bilities in the acceleration voltages, resulting in a spread 
of individual ion energies Ei. We approximate the ion 
energy distribution Pib(ifi, Ajb) by a symmetrical energy 
spread function peaking at Ajb with a FWHM of Ajb. 

The voltage for the first-stage acceleration on the high- 
voltage platform was stable to within Ri 5 V. For the 
second-stage acceleration by the high-voltage platform 
potential, two different voltage sources have been used. 
The first voltage source, hereafter denoted as A, was lim¬ 
ited to 40 kV and was used for most measurements with 
Aib < 53 keV. This source was stable to within r: 4 V. 
The spreads from the first- and second-stage accelera¬ 
tion, added in quadrature, result in a total ion energy 
spread of r 6 eV. As will be shown in Section PV A1 
this ion energy spread is negligibly small in comparison 
to the detector resolution. 

For Aib > 53 keV a less stable voltage source B was 
employed in the second-stage acceleration. The voltage 
fluctuations here dominate over those from the first-stage 
acceleration. Direct determination of the voltage spread 
was difficult and therefore we derived it from the differ¬ 
ence in peak widths in the detected energy spectra 
and P® acquired with the same ion type, the same Aib, 
but with voltage sources A and B, respectively. As shown 
below in Section IIV A1 the spectra for atomic projectiles 
are dominated by a single Gaussian-like peak. Thus, by 
approximating Pib(£’i, Pib) with a normal distribution we 
estimated the ion beam energy spread for the voltage 
source B as 



where A^ and A® are the FWHMs of P^ and P®, re¬ 
spectively. We applied the procedure separately for the 
H+ and Ar+ data acquired at ^ 53 keV and obtained 
values of A® which are equal to within their statisti¬ 
cal uncertainties. The error-weighted average of these 
two results is A^ = 320 eV. A conservative estimate of 
the total systematic error of this procedure is ±100 eV, 
which exceeds the statistical error by a factor of ^ 2. In 
Section [IV Al we discuss how we are able to improve our 
estimate for the systematic error. 


B. Detector setup 

The detector setup was developed at the Kirchhoff In¬ 
stitute for Physics of the University of Heidelberg, Ger¬ 
many. The MMC detector itself was derived from the 
well-established maXs-200 design, which was originally 
developed for X-ray spectroscopy^. For the present mea¬ 
surements we have used only two detector pixels, each 
having a 0.45 x 2.0 mm^ sensitive area. The absorber for 
one of the pixels was made of a 50 ^m thick gold layer. 
If not explicitly mentioned otherwise, the data presented 
below have been obtained from this pixel. Selected mea¬ 
surements were performed with the second pixel where 
the 50 /rm thick gold absorber was coated with an ad¬ 
ditional 180 nm thick layer of aluminum. In both cases 
incident atoms and molecules were stopped in the upper¬ 
most ^150 nm layer of material, as verified by calcula¬ 
tions using the SRIM software (“Stopping and Range of 
Ions in Matter”—). Hence, the two detector pixels corre¬ 
spond to stopping in gold and in aluminum, respectively. 
The absorber itself has full contact to a 2.4 /rm thick 
thermal sensor made of a dilute gold-erbium alloy, mixed 
at a mole fraction of [Er]/[Au±Er]=885 ppm. The er¬ 
bium was depleted in ^®^Er, which has a large hyperfine 
splitting and would otherwise add a large heat capacity 
at the operational temperature. A detailed description 
of the maXs-200 microfabrication and performance has 
been given previously^2i^. 

The ~ 10 mK working temperature of our MMG detec¬ 
tor was maintained with a dilution refrigerator. To min¬ 
imize the heat flux from the 300 K black-body radiation 
of the ion beamline onto the detector, a set of radiation 
shields at 800 mK, 4 K, and 50 K surrounded the detec¬ 
tor. Traveling along the ion beam trajectory, the beam 
encountered a set of sequential ~ 1 cm apertures at 50 K 
followed by a single 4 K foil with ~ 4 /rm-diameter pin¬ 
holes to limit the radiation flux onto the detector head 
to below 1 nW. The pinholes were located so that there 
was at most one in front of each detector pixel. The pixel 
used in the analysis was selected by connecting the read¬ 
out electronics described below. The rate of detected 
events was maximized before each data run by slightly 
steering the ion beam. 


C. Data acquisition and processing 

For each combination of ion species and ion beam en¬ 
ergy we accumulated between 500 and 2000 detection 
events. For each event the electronic pulse from the ther¬ 
mal sensor was read out via the SQUID electronics, am¬ 
plified, and digitized. The individual pulses were ana¬ 
lyzed to derive the detected energies and thus also 
the energy distribution P^{E^, Ajb, M). 

The electronic signal of the detector reflects the tem¬ 
perature of the magnetic moments of the paramagnetic 
sensor, which closely follows the absorber temperature. 
The temporal evolution of the signal is determined by 
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the thermal links between the absorber, the sensor, and 
the thermal bath. In our setup a typical pulse reached its 
maximum amplitude within ~ 30 /rs after the detection. 
The subsequent cooling results in an exponential signal 
decay with a time constant of ^ 5 ms. In our experiments 
the energies detected were small enough so that the sen¬ 
sor heat capacity can be approximated by a constant over 
the whole range of sensor temperatures occurring during 
the detection event. Additionally, the detector-operation 
temperature range was chosen so that the sensor magne¬ 
tization was proportional to the sensor temperature. As 
a result, the shapes of the acquired electronic pulses were 
independent of the detected energy. Moreover, to a very 
good approximation the pulse amplitudes were linearly 
proportional to the heat released in the absorber. To ob¬ 
tain the detected energies on a relative scale we fitted the 
pulse amplitudes using a fixed pulse shape. Specifically 
we used the advanced fitting procedure based on optimal 
filtering, also called “matched filter”— 

To put the detected kinetic energies on an absolute 
scale we additionally detected photons from an ^"'’^Am 7 - 
source placed at ^ 15 cm distance from the detector. 
The photon energy from such a source peaks sharply at 
59.5409 keV— . In our setup we have determined the 7 
energies with ^15 eV precision, which corresponds to 
^ 0.25% uncertainty in the scale factor of the energy 
calibration. The non-linearity of the energy scale is ex¬ 
pected to be less than 0.1% at 60 keV— . 

D. Detector energy-response function 

The detector energy-response function relates the im¬ 
pacting ion kinetic energy Ei to the measured energy . 
As explained in Section III B[ the energy-response func¬ 
tion of an MMC detector setup is determined by statis¬ 
tical processes. First, thermal processes in the absorber 
transform the kinetic energy Ei into heat Ah with a sta¬ 
tistical distribution Fj^{Eh, Ei, M). This process may 
depend on the type of incident particles denoted by M. 
Second, the thermal noise in the detector and the vari¬ 
ous noises in the readout system convert Ah into the de¬ 
tected energy Am. The response function corresponding 
to this second transform, A°gj(Am,Ah), reflects the in¬ 
trinsic detector resolution. The overall detector response 
function is then a convolution of the two spreads in the 
Ah-domain: 

Adet(F;m,i^i,M) = A°et(Am,Ah)®A7P(Ah,Ai,M). (5) 

We show below that the detector energy-response func¬ 
tion Adet can have a complex shape, but is still often 
dominated by only a single peak. When that is the case 
we can define the detector resolution Adet(^'i, M) as the 
FWHM of Adet(^^m) for a given Ei and projectile M. The 
relative resolution is then defined as Adet{Ei,M)/Ei. 

We have determined the intrinsic-resolution function 
of our detector setup, A°gj, by applying our standard 
event analysis procedure (as described above) on the 


time-dependent electronic signal acquired between the 
detection events. Fluctuations of this baseline signal arise 
from the same instabilities and noise sources in the detec¬ 
tor and the data acquisition system that also act on the 
particle-generated pulses. As explained above, the detec¬ 
tor signal scales linearly with the energy input. Thus the 
pulses derived from this baseline fluctuation reflect the 
intrinsic noise of the detection system and the distribu¬ 
tion of their peak amplitudes can be used as a measure 
for the intrinsic energy resolution of our detector setup. 
The inferred A°gj displayed a Gaussian shape centered at 
0 eV with a FWHM of = 107 ± 5 eV. 

IV. RESULTS AND DISCUSSION 

We have acquired MMC energy spectra Am(7?m) for 
various singly charged atomic and molecular ions in the 
energy range from Aib ~ 13 to 150 keV. We first present 
our measurements with atomic ions in Section IIV Al In 
Section llVBI we demonstrate the mass spectrometric ca¬ 
pabilities of the MMC detector using molecular fragmen¬ 
tation measurements. In Section IIV Cl we discuss the 
origins for the various features in the acquired spectra. 
The data in these first three subsections were collected 
using the uncoated gold absorber. In Section IIV PI we 
present the effects of the MMC absorber material on the 
shapes of the detected spectra. Then, in Section ITV El we 
present exploratory observations revealing neutral frag¬ 
mentation products from the acetone radical cation, the 
largest molecule studied here. 


A. MMC detection of atomic ions 

In Figure [1] we present four 7?ib,H+) spectra 

acquired with proton beams at Aib = 14.7, 52.6, 89.9, and 
151.5 keV. In the figure we label the spectra as Ha, Hb, 
He, and Hd, respectively. Each spectrum is dominated 
by a sharp peak close to Aib. A fit of the spectra in 
the vicinity of the peak using a Gaussian function gives 
EWHMs of Am = 117 ±7, 165 ±29, 401 ±10, and 415 ± 
26 eV. The errors given are the statistical uncertainties 
only. 

Each of the acquired spectra is a convolution of the 
detector energy-response function and the ion energy dis¬ 
tribution: 

±m(Ai„, Aib, fW) = Fdet(75mj Ai, Af) 0 Pib(Ai, Aib). (6) 

Eor the Ha and Hb spectra, voltage source A has been 
employed to accelerate the proton beams. In this case, 
the ion energy distribution width A(^, which has been 
determined independently (see Section IHI Al) . is negli¬ 
gibly small compared to the intrinsic detector resolution 
A°et derived in Section lTHDl Hence, the ion beam energy 
distribution Pib can be neglected in Equation dH). Cor¬ 
respondingly, the Ha and Hb spectra reflect directly the 
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FIG. 1. The MMC-detected energy spectra for proton beams at = 14.7, 52.6, 89.9, and 151.5 keV are plotted in gray as 
a function of the detected energy Em- We label the spectra as Ha, Hb, He, and Ha, respectively. The values of Pm are given on 
an arbitrary scale. In order to visually separate the spectra, vertical offsets have been added to Ha, Hb, and He. The 7 -source 
signal is labeled in He. The insets display the spectral details close to the respective beam energies Pib. The statistical error 
bars are displayed only in the insets. Gaussian-fit curves of the dominant peak in each spectrum are plotted in blue. The blue 
horizontal solid lines indicate the respective FWHMs and the dashed lines the baselines. 


overall energy-response function Pjet and the correspond¬ 
ing overall detector resolution is given by Adet = Am. 
On the other hand for the spectra He and Hd, acquired 
with voltage source B, the widths of the prominent peaks 
are larger by a sizable factor, in spite of the only slightly 
higher ion energies. We conclude that the widths of these 
spectra are strongly influenced by the ion beam energy 
spread. From separate measurements (explained in Sec¬ 
tion |llT^ we derived a value of Af^ = 320 ± 100 eV for 
this energy spread. From the measured widths, and ap¬ 
proximating both the measured spectra Pm(^'m) and the 
ion energy distribution Pib(£'i) as normal distributions, 
we obtain the detector resolution as 

Adet « V^(Am)2-(AB)2. (7) 

For the He and Hd spectra, we find Adet = 241 ± 17 eV 
and 264 ± 41 eV, respectively. The given errors reflect 
only the statistical uncertainties of the Am widths. The 
inferred width Adet resulting from Equation Q must 
be larger than the intrinsic detector resolution A^^^ = 
107 eV. This requirement and the measured Am for spec¬ 
tra He and Hd set an additional limit on the ion energy 
spread to be Af^ < 386 eV. This enables us to reduce 
the estimated systematic uncertainty on the beam en¬ 
ergy spread A|^ to t igg eV. 

In addition to protons, we have studied the MMC reso¬ 
lution for the atomic ions C’*', N+, and Ar+. In Figure [2] 
we compare the corresponding energy spectra Pm(Em) 
acquired at beam energies around Eib = 53 keV. Similar 
to the proton spectra, the spectra of the heavier atomic 
ions are also dominated by a peak close to Eib. In con¬ 
trast, however, this main peak now has a clear tail ex¬ 
tending towards low energies, suggesting that a small but 


significant part of the ion kinetic energy was not detected 
for a noticeable fraction of events with heavier atomic 
ions. In order to quantitatively describe the resulting 
asymmetric peak shape we fit the spectra by a convolu¬ 
tion of a Gaussian function and a left-sided exponential 
function^iH 


/(^m) =^exp 
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where A is the peak area, t is the exponential decay pa¬ 
rameter of the left-sided asymmetric part, Eq is the peak 
position, a is the Gaussian width, and erfc is the comple¬ 
mentary error function. The fitting range was iteratively 
adjusted to Eq — 3a — 3t < < Eq + 3a. 

In Figure [3] we present the peak widths A^et and the 
asymmetry parameters r/A^et for all of the atomic ion 
spectra. The measured width of each spectrum, A„i, is 
the numerically determined FWHM of the fitted function 
f{Eni). As explained above for protons, the detector res¬ 
olution for measurements with voltage source A is taken 
as Adet = Am. Similarly, we have corrected the widths 
Am for spectra acquired with the voltage source B using 
the approximate method given by Equation 0. 

The derived detector resolutions Adet for atomic ions. 
Figure [31(a), display a clear dependency on the ion type. 
The quantity Adet grows with increasing ion mass. The 
best MMG energy resolution was achieved for proton 
beams with Adet < 200 eV at all measured values 
of Ejb. The narrowest spectrum was found with pro¬ 
tons at Aib = 14.7 keV, giving Adet = 117 ± 7 eV. 
The best relative resolution was achieved with protons 
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FIG. 2. The MMC-detected energy spectra Pm for H+, C^, 
N^, and Ar'*' beams at beam energies aronnd Pib = 53 keV 
are shown by gray data points in panels (a) to (d), respec¬ 
tively. The data for H’^ and Ar"*' were acquired using voltage 
source A, while data for C'*' and were acquired using volt¬ 
age source B. The statistical uncertainties are given by the 
vertical error bars. Green arrows mark the respective nomi¬ 
nal ion beam energies Pib. Fits by Equation (ISJ are shown 
by the blue curves. The blue horizontal lines indicate the 
respective FWHMs. 


at Pib = 151.5 keV, namely Auet/T^ib ~ 1-7 x 10“^. 
The peak widths increase for heavier atoms, reaching 
Adet ~ 700 eV for Ar+. Given the large statistical uncer¬ 
tainties and the systematic uncertainty in A^ (relevant 
for the hollow symbols), no clear dependency of Adet on 
Pib can be identified. 

In Figure [3Kb) we show the peak asymmetry described 
by the ratio of the left-sided exponential decay constant 
T to the peak FWHM (i.e., Adet)- The most symmetric 
spectra - nearly Gaussian - were acquired with protons. 
Fitting these spectra using Equation ([ 8 ]) failed in some 
cases due to numerical divergence at r —> 0. In such cases 
we have used a pure Gaussian to fit the data and report 
a generic value for r of “0”. Heavier atoms, compared 
to the protons, display more asymmetric peaks with a 
slightly decreasing tendency for increasing P^. 

In all of the atomic ion measurements, most events 
were detected within the main peak, i.e., at energies of 
Pq — 3cr — 3 t < Pm < Po + 3cr. At higher energies no 


counts were recorded, which demonstrates the extremely 
low-background detection capability of MMC detectors. 
In the low-energy tails, the spectral count rates are small, 
decreasing smoothly towards lower energies (ignoring the 
events at ^ 59.5 keV from the 7 source). In Figure |4| we 
plot the fractions of event counts detected in these far 
low-energy tails, defined as 

‘ J,°^Pm{Em)dE^ ■ 

Rt reaches values as low as ~ 10% for high-energy proton 
projectiles. This fraction grows to up to ^ 35% for heav¬ 
ier atoms. For H+ and G'*' a clear decrease of Rt with 
increasing Pjb can be seen. This trend cannot be seen in 
the N+ data and even seems to be inverted for Ar+ ions. 
Measurements presented in Section IIVDI show that the 
tail fraction is significantly reduced with the Al-coated 
absorber. 


B. MMC mass spectrometry on molecular ions 

In addition to atomic ions we have also studied the 
calorimetric detection for several types of molecular ions. 
As an example we plot the detected energy spectrum for 
a beam of GH;^ at Pjb = 150.6 keV in Figure [5] (black 
line). Similar what we found for atoms, this molecular 
spectrum is dominated by a single peak near to Pib. From 
the peak width of Am = 786 ± 71 eV and from the ion 
beam energy spread A^®, we derive a FWHM detector 
resolution of Adet = 718 ±78 eV. This exceeds the energy 
width for G"*" atomic ions (Figure [3]) by ^ 200 eV. 

Using an energetic molecular ion beam allowed us to 
demonstrate the MMG mass spectrometry capabilities di¬ 
rectly in a single measurement. For this we directed all 
of the charged particles away from the detector, using 
a magnetic deflector placed ^ 1 m upstream of the de¬ 
tector. In this configuration, the only remaining particles 
hitting the detector were neutral atoms and molecules re¬ 
sulting from collisions of the molecular ions with residual 
gas in the ion beam line. Clearly, only a small fraction of 
the collision-induced neutrals are able to reach the MMC 
detector through the narrow pinhole in the mask in front 
of the detector. Thus, to reach acceptable data quality 
from this measurement configuration we have used high 
ion beam intensities and increased residual gas pressure 
in the beamline. Processes leading to neutral fragments 
are discussed in Section IIV El 

The red curve in Figure jS] displays the energy spectrum 
of neutral fragments resulting from collisions of CH;^ ions 
at Pib = 150.6 keV with residual gas molecules along 
the upstream beam line and reaching the detector pixel 
through the pinhole mask. The spectrum displays dis¬ 
tinct peaks at energies matching well to the respective 
fragment-to-parent mass fractions of Pib for all possible 
fragments of the CH;^. The peak widths in the neutral 
spectrum are A^ = 333 ±20 eV for H, 349 ±46 eV for H 2 , 

























FIG. 3. Panel (a): The widths Adet of the dominant peak feature in the MMC-detected energy spectra for various atomic 
ions are displayed as a function of the ion beam energy Eih- The legend for ion types is displayed on the very right-hand 
side. The ion beams for the full-symbol data were accelerated by voltage source A. The correspondingly low ion energy spread 
was neglected. The hollow symbol data were acquired with a less stable beams using voltage source B. These widths were 
corrected for the ion energy spread using Equation 0. The error bars reflect the statistical uncertainties of the fitted widths, 
but do not include the eV systematic error on the ion energy spread. The dashed line indicates the average intrinsic 
detector resolution of = 107 eV. Panel (b): The peak asymmetry expressed as the ratio of the left-sided exponential decay 
parameter and the peak width (r/Adet) is plotted as a function of Eih- The meaning of the full/hollow symbols and of the 
error bars is identical to panel (a). The zero amplitude of some points (protons only) indicate fully symmetric peaks where the 
fit by Equation failed for numerical reasons. 



FIG. 4. The fraction of events detected at energies below 
the dominant peak, Rt, is plotted as a function of ion beam 
energy Eih for various atomic ion beams. The exact definition 
of the Rt fraction is given in the text. The error bars indicate 
statistical uncertainties. 


908 ± 180 eV for C, 1230 ± 67 eV for CH, 1013 ± 155 eV 
for CH 2 , and 813 ± 78 eV for CH 3 . The CH 3 peak width 
is equal to that from the CH^ spectrum (to within the 
statistical uncertainty). This suggests that the molecular 
velocity is not significantly altered in the non-dissociative 
neutralizing collision with the residual gas. Moreover, 
this also allows us to conclude that the charge of the 
detected particle is irrelevant for the calorimetric mea¬ 


surement, to within the detector resolution. 


Interestingly, the peaks for the CH and CH 2 fragments 
are clearly broader than that for CH 3 . The additional 
spread may originate from the kinetic energy released in 
the CH 3 dissociation process. For example, for a 150 keV 
CH 3 beam, a center-of-mass kinetic energy of 1 eV re¬ 
leased into CH and H 2 fragments dissociating along the 
ion beam axis transforms into an ~ 260 eV change of the 
fragment kinetic energy in the laboratory frame. Since 
the present measurements detect only a small fraction 
of the neutral fragments produced in the beam, the rel¬ 
ative detection yields for the different fragment species 
may depend significantly on the energy released and the 
angular emission pattern of the respective fragmenta¬ 
tion channels. Nevertheless, we find that even under 
these non-ideal conditions the high-resolution calorimet¬ 
ric mass spectra still yield qualitative information about 
the molecular composition and the fragmentation chan¬ 
nels. 


In a more pronounced fashion than for the atoms, the 
mass peaks for both the primary molecular ion and the 
neutral molecular fragments are accompanied by tails 
that will be further discussed below. These wider tails 
turn out to be significantly influenced by the coating of 
the absorber in which the molecular stopping occurs (see 
Section HV DU . 
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FIG. 5. The MMC-detected energy spectra for a CH^ beam at i?ib = 150.6 keV is plotted by the black line as a function 
of the detected energy E^. The red curve plots the spectrum acquired with ions prevented from reaching the detector, so that 
only neutral CH 3 and its neutral fragments were detected. The data have been vertically shifted from zero in order to visually 
separate the two spectra. The rightmost green arrow indicates the nominal ion beam energy i?ib and the remaining green 
arrows, going from left to right, show the energy fraction corresponding to the masses of H, H 2 , C, CH, and CH 2 fragments. 


C. Modeling MMC energy spectra 

As was shown above fSection llV Ajl . the detected en¬ 
ergy spectra for atomic projectiles are dominated by a 
single narrow peak close to the nominal ion beam energy 
Aib, while more complex shapes are found for molecu¬ 
lar projectiles. Thus, the CHg data at £lib = 150.6 keV 
(FigurelU black line) display a low-energy shoulder on the 
main peak in the energy range from ^ 140 to 147 keV. We 
find that the relative intensity of this feature grows sig¬ 
nificantly for lower ion beam energies, as seen in the CH^ 
data acquired at Aib = 53.3 keV (FigurelHl). Similarly, we 
also observed prominent low energy shoulders in energy 
spectra recorded for CO^ and (CH 3 ) 2 CO+ beams. 

In order to understand the detected spectral shapes we 
have modeled the processes occurring during the stopping 
of molecular projectiles in the MMC absorber. In our 
simplified model we assume that the detected molecule 
breaks into separate atoms immediately after impinging 
on the detector surface and neglect the few eV of molecu¬ 
lar binding energy. Additionally, we disregard the charge 
state of the projectile particles, based on our findings de¬ 
scribed earlier. With these assumptions we can study 
stopping processes in the absorber separately for each 
of the atomic constituents of a given molecule, such as C 
and H in case of CH;^. The model for the molecular spec¬ 
tra can then be constructed from a suitable combination 
of the modeled atomic spectra. 

We have employed the SRIM software to simulate the 
energy loss processes of the projectile atoms in the MMC 
absorber. SRIM^ is a 3D Monte-Carlo code that simu¬ 
lates the propagation of ions in matter. It describes not 
only collisions between the projectile and target atoms, 
but also accounts for target-target collisions due to cas¬ 


cades triggered by recoiling target atoms, backscattering 
of projectiles, and sputtering of atoms from the target 
material. The resulting detailed data on each collision 
allows us to evaluate the amount of projectile kinetic en¬ 
ergy which is not transformed into heat and thus escapes 
calorimetric detection. Specifically, for each simulated 
detection event we have calculated whether the projec¬ 
tile atom was back-scattered, and if yes, then what frac¬ 
tion of the kinetic energy has remained in the absorber. 
Next, we have evaluated the number and the total ki¬ 
netic energy of target atoms sputtered from the detector 
surface. Lastly, for each impact we track the number of 
defects created in the absorber lattice. Such defects arise 
when an atom of the target lattice is kicked out by re¬ 
coil from its stable position (either by the projectile atom 
or another recoiled target atom). The resulting pair of 
an interstitial atom and a vacancy in the lattice (i.e., a 
Frenkel pair or a Frenkel defect^) holds a potential en¬ 
ergy i^FP. This non-thermal energy storage results in 
a reduction of the energy detectable by the MMC. The 
SRIM code calculates the number of Frenkel defects along 
each incident ion track. To obtain the total energy loss 
due to Frenkel pair formation we multiply the number of 
Frenkel pairs by App which is a parameter of our model. 
The number of replacements where an energetic target 
atom terminates its movement in a previously created 
vacancy, is negligible. Altogether, for each impact event 
we have subtracted the above listed energy losses from 
the initial kinetic energy of the incident atom and thus 
obtain individual atomic detected energies E^°. For a 
given molecule, we simulate 10^ impact events related 
to each atomic constituent type. This yields the model 
energy spectra for each atomic species. 

An additional loss of energy may arise from secondary 
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FIG. 6. The MMC-detected energy spectra for CH^ beam at i5ib = 53.3 keV is plotted by the black line as a function of 
detected energy j5m. The red curve represents our SRIM-based model spectrum Pm, assuming a Frenkel-pair energy Pfp = 
2.8 eV. Going from right to left, the green vertical lines mark the nominal ion beam energy Pib and the maximal energy losses 
due to backscattering of one, two, and three H-atoms, and one C-atom. The model spectrum is scaled in amplitude so that the 
integral between Pm = 5 keV and 54 keV matches the experimental data. 


electron emission by the detector absorber after the im¬ 
pact of a particle. Studies on secondary electron emis¬ 
sion from solid gold targets, however, have found only 
small secondary electron yields of < 2 electrons per in¬ 
cident ion^iS and mean secondary electron energies of 
^10 eV— for a broad range of projectile energies and 
projectile types. Based on those studies we estimate the 
energy loss due to secondary electron emission to be neg¬ 
ligible compared to the other energy loss processes and 
we do not include it in our model. 

To calculate the detectable energy for a molecular pro¬ 
jectile, P“°*, we have randomly picked simulated events 
for each of the constituent atoms and sum the corre¬ 
sponding E^° values. Then we have added a random 
Gaussian spread describing the intrinsic detector res¬ 
olution = 107 eV) and the ion energy spread 

(A^ = 320 eV FWHM), obtaining the model energy dis¬ 
tribution PmiEm)- Adding the ion beam energy spread 
only in the last step (instead of varying Ei as the input 
parameter of the SRIM simulations) is justified by the 
small relative ion beam energy spread of A^/Pib < 1%, 
within which the relative contributions of the various en¬ 
ergy loss processes stay nearly constant. 

In Figures|n]and[3we compare the model energy distri¬ 
bution Pyn for CHg at Pib = 53.3 keV to the correspond¬ 
ing experimental data. In red we plot the result of the 
complete model including all of the energy loss processes 
mentioned above and for the value of Ppp yielding the 
best agreement with the data. Partial effects are shown 
in Figure[3for the Frenkel-pairs only (blue) as well as for 
Frenkel pairs plus sputtering (green). Comparing this 
latter curve with the full result (red) shows the rather 
strong effect of the backscattering process. The curves 
for the Frenkel-defects only are also shown for a range of 


Ppp values. 

The shape of the dominant peak at E^ = 52.7±0.5 keV 
can be assigned to the energy spread due to Frenkel pair 
formation. The Frenkel-pair energy used in the model 
affects only this feature in its width and mean energy. 
Comparing directly the peak energy in the model to the 
measured data is difficult due to the systematic uncer¬ 
tainty in the nominal ion beam energy Pjb (Sec. IIII Al) . 
Instead we have varied the Ppp value in steps of 0.05 eV 
and in each iteration compared the width of the domi¬ 
nant peak in the model spectra to that in the measured 
data. Simultaneously we adjusted the model beam en¬ 
ergy so that the center of the peaks overlap. The best 
match was reached for Epp = 2.8 eV and by adjust¬ 
ing the model beam energy by -1-0.1 kV (this is only 
0 .2% of Aib) with respect to the nominal beam energy 
of Aib = 53.3 keV. Such an energy shift is well below 
the estimated Aib uncertainty. The energy losses due to 
sputtering of absorber atoms add only a small left-handed 
tail on the main peak and do not produce any distinct 
spectral feature. Backscattering of H atoms is responsi¬ 
ble for the wider feature between ~ 45 and 52 keV and 
backscattering of C for the tail remaining visible down 
to ^ 15 keV. The small feature at ^ 1.2 keV was not 
observed in most of other acquired spectra and is inter¬ 
preted as a spurious event not originating from CHj^. 

Backscattering appears to be the key for explaining 
the differences between the atomic and molecular spec¬ 
tral shapes beyond the main peak. For atomic projectiles, 
the observed low-amplitude and low-energy tails are at¬ 
tributed to backscattering of the projectile atoms them¬ 
selves; any of the back-scattered atoms can take away up 
to the full original kinetic energy. Typical spectra can 
be seen in Figure [TJ On the other hand, for molecular 
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FIG. 7. A comparison of the MMC detected energy spectrum (gray) and a corresponding SRIM-based model for CH^ (red). 
The thick blue line plots the model spectrum assuming energy losses only via Frenkel-pair formation with an adjusted Frenkel- 
pair energy of Eyp = 2.8 eV. The thin blue lines plot the same model using an alternative value for F'fp of 6.0 eV, 4.0 eV, and 
2.0 eV, respectively. The green line plots the model assuming Frenkel-pair formation and sputtering. The difference between 
red and green line then reflects the effect of backscattering. The horizontal green bar labels the uncertainty range for the 
nominal ion beam energy i?ib = 53.3 ± 0.4 keV in the experiment. The red arrow marks the adjusted ion beam energy for the 
model. 


projectiles the dissociated atomic fragments are expected 
to backscatter independently, each fragment being able 
to carry away only the part of the original projectile ki¬ 
netic energy given by the fragment-to-molecule mass ra¬ 
tio. In the case of CH^, energy loss maxima of Ew^/lb 
and of 12iiiib/15 result for backscattered single H and C 
atoms, respectively. Given the H multiplicity, the most 
likely backscattering process is that of a single H frag¬ 
ment, taking away energy between 0 and ^ iiiib/15. As 
a consequence we observe a strong, broad feature in Pm 
between Pm ~ and Pm ~ Pib- Structures 

in Pm at lower energies decrease in magnitude as the si¬ 
multaneous backscattering of multiple H atoms becomes 
less probable. The lowest energies can be accessed only 
by backscattering of carbon atoms. 

Only a few parameters enter our model. The most 
influential one is the Frenkel-pair energy which directly 
scales the width and position of the main peak. Our 
fitted value Ppp =2.8 eV has a strong correlation with 
the beam energy spread A®. The tioo eV uncertainty 
of A® causes an uncertainty of Iq ® eV in the adjusted 
value of Ppp. 

Another important parameter entering the SRIM 
model is the threshold displacement energy Pdispb i-e., 
the minimum energy needed for a recoiling target atom 
to overcome the lattice barrier, move away from its orig¬ 
inal position, and form an interstitial. In our model we 
use Pdispi = 43 eV, a value recommended in a recent 
review^. The Pdispi values listed in that review range 
from 30 eV to 44 eV. Using this range as the uncertainty 
in Pdispi results in an uncertainty of (to ’7 eV for Ppp. 

Furthermore we set the lattice binding energy to 
^^latt = 3.81 eV— and the surface binding energy for 


gold to Psurf = 3.83 eV— . In the SRIM simulation Piatt 
is the energy that every recoiling target atoms loses when 
it leaves its lattice site. Psurf is relevant for calculating 
sputtering. The estimated 0.5 eV uncertainties on Piatt 
and Psurf propagate to a negligible uncertainty in Ppp. 
For details on the relations between Pdispi, i'latt, and 
Ppp see, e.g., the SRIM documentation^ or Figure 11 
m^. 

Adding all the systematic uncertainties together in 
quadrature gives a Frenkel-pair energy of Ppp = 2.8 ± 
0.9 eV. To the best of our knowledge, no other experi¬ 
mental values of Ppp have been published for gold. Cal¬ 
culated values range from Ppp = 3.0 to 4.1 eV, depend¬ 
ing on the calculational method as well as on the specific 
lattice defect type formed^iiZ^. The Frenkel-pair energy 
determined from our model matches well with these cal¬ 
culated values. 


D. Optimizing MMC resolution - absorber materials 

Practical usage of the MMC detector for spectrom¬ 
etry requires a narrow energy-response function. We 
have shown that the dominant peak in all of the col¬ 
lected energy spectra is narrow and that relative reso¬ 
lutions down to ^ 10“^ can be reached. Nevertheless, 
the low-energy tails of the energy-response function can 
significantly complicate the analysis of energy spectra, 
especially for molecular projectiles. 

Our modeling studies indicates that the low-energy 
tails in the measured spectra originate predominantly 
from backscattering. Thus, in order to optimize the 
MMC energy-response function, we have searched for 
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FIG. 8. MMC detected energy spectra for CH^ at Eih = 
53.3 keV, acquired with the gold-absorber pixel (black) and 
with the aluminum-coated pixel (red). The amplitude of the 
spectra have been scaled such that the total integrals are 
equal. An artificial offset in Pni has been added for clarity. 
The green arrows label the nominal ion beam energy Eih- 


absorber materials with low backscattering probabilities. 
SRIM models show a general trend of lower backscatter¬ 
ing for materials composed of lighter atoms (i.e., lower 
atomic number Z). Backscattering is determined pre¬ 
dominantly by the first few surface monolayers. Hence 
only a thin absorber coating of a low-Z material should 
be sufficient to improve the MMC energy-response func¬ 
tion. Moreover, SRIM calculations indicate that fewer 
Frenkel pairs should be created in \ow-Z materials. Thus, 
provided that the Frenkel-pair energy does not depend 
on the absorber material or is lower for a specific \ow-Z 
material, then the dominant peak should also narrow. 

To test this hypothesis we have coated the 50 ^m 
gold absorber of one of the MMC pixels with 180 nm 
of aluminum. In Figure [5] we compare the energy spec¬ 
tra from the bare and Al-coated gold pixels, acquired 
with the same beam of CH;^ at iJib = 53.3 keV. Clearly, 
using an A1 coating results in a strong suppression of 
the backscattering. There is a clear reduction near 50 
keV in the shoulder due to backscattering of H atoms, 
and also a clear decrease below 46 keV in the tail due 
to backscattering of C atoms. Interestingly, the main 
peak is shifted towards lower energies when using the 
Al-coated absorber. This cannot be reproduced by our 
SRIM-based model. We anticipate that the shift is due to 
the aluminum being superconductive at the MMC oper¬ 
ating temperatures, which is below the aluminum critical 
temperature of Tc = 1.2 K. Hence, the energy loss pro¬ 
cesses in the Al-coated pixel may be different from those 
for normal-conducting gold. This interpretation is sup¬ 
ported by previously published experiments investigating 
the properties of various superconductors when used as 


a particle or X-ray absorber in a micro-calorimeter. At 
first glance, superconductors seem to be a very promis¬ 
ing class of cryodetector materials, considering their rel¬ 
atively small specific heat at temperatures well below T^. 
However, as shown byS, which summarized earlier work 
on various superconductors, and by more recent stud- 
ies2i“— on aluminum and rhenium, the thermalization of 
energy in superconductors is not well understood at tem¬ 
peratures as low as that of our absorber at ~ 10 mK. At 
these temperatures there is very little phonon excitation 
(T < 2 X 10“"^ 0D, where 0d is the Debye temperature) 
and surprisingly long thermal pulse decay times have 
been observed^^. Analysis of these measurements sug¬ 
gests that a significant fraction of the deposited energy, 
from a few percent to tens of percent, could be stored in 
very long-lived excitations with lifetimes ranging up to 
several seconds. This fraction of energy would escape de¬ 
tection. The relevant long-lived excitations could be due 
to quasi-particles resulting from broken Cooper pairs, ex¬ 
citations of magnetic flux lines, or excitations of nuclear 
spins or atomic tunneling systems, to name just a few 
possibilities. In fact, the apparent magnitude of the en¬ 
ergy loss in our detector with the aluminum coating, on 
the order of 5%, fits well to these previous observations 
and interpretations. 

We also used the Al-coated detector pixel to reproduce 
the measurements on neutral CH;^ fragments created by 
residual gas collisions of CHg ions at i?ib = 150.6 keV 
(see Section ITVBI) . The results are displayed in Figure IH] 
A comparison to the data acquired with the bare gold 
absorber (Figure [5]) shows, again, a strong reduction of 
the low-energy tail when the Al-coated absorber is used. 


E. Acetone radical fragmentation 

Using our MMC setup we have also investigated the 
fragmentation of acetone radical cations, (CH 3 ) 2 CO+ 
(58 u), due to residual gas collisions. To detect the neu¬ 
tral fragments we followed the approach outlined in Sec¬ 
tion |TVBl The data acquired at Uib = 150.45 keV using 
the Al-coated detector pixel are displayed in Figure [TUI 
Similar to our other measurements with this pixel, the 
measured energies are ^ 4% lower than those obtained 
with the gold-only pixel. For this large molecule, the data 
obtained with the uncoated gold pixel (not shown here) 
are much less resolved than those from the Al-coated 
pixel, although, as was consistently seen, the full-mass 
peak for the uncoated pixel lies much closer to Ujb. In 
the neutral fragment spectrum of Figure (TU] (collected 
using the Al-coated pixel) we assigned the most distinct 
peaks by multiplying the fragment-to-molecule mass ra¬ 
tios by a full kinetic energy of only 96% of Uib (the en¬ 
ergy step corresponding to 1 u is then 2.49 keV). The 
small shifts between the peak centers and the predicted 
fragment energies, seen for only a few specific fragment 
masses (notably near 13 u at ~ 32 keV), may be due to 
more complex, still unexplained effects in the stopping of 
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FIG. 9. Same as Figure [S] but acquired with the Al-coated MMC absorber. 


the molecular species. 

Two processes are expected to generate the neutral 
fragments. The first is collision induced dissociation 
(CID, also known as collision activated dissociation - 
CAD), which does not change the charge of the target 
particle (a residual gas molecule in our case). In CID 
at least one of the projectile fragments stays charged. 
The second process is an electron transfer from the 
neutral target to the projectile cation neutralizing it. 
As in electron-ion recombination, the ionization energy 
released then usually leads to the dissociation of the 
neutralized projectile molecule into neutral fragments 
(known as electron transfer dissociation - ETD). The 
electron transfer may also result in an internally excited, 
but non-dissociated, neutral molecule. 

We also want to emphasize here that the setup used in 
these exploratory studies was not designed for quantita¬ 
tive branching ratio measurements. The kinetic energy 
released in the various CID and ETD outgoing channels 
may vary significantly. This affects the transverse distri¬ 
bution of the fragment impact positions on the detector. 
As a result, the true geometrical efficiency for detecting 
the fragments through the pinhole in front of the detec¬ 
tor may vary strongly as well. Thus, the relative peak 
heights should be interpreted with care. Nevertheless, we 
find that the high mass resolution of our setup provides 
unprecedented insight in the CID and ETD channels pro¬ 
ducing neutral fragments by showing a number of discrete 
fragment mass peaks. 

Many studies exist on CID and ETD of the acetone 
radical2i“— . However, the vast majority of the tech¬ 
niques used collected only charged reaction products. 
Given this limitation, these studies assigned the domi¬ 
nant dissociation channel to 

(CH3)2CO+-bA^CH3CO+-bCH3-bA. (10) 

Here A is the target particle, which in our case of a non- 
baked UHV beam line is expected to be H 2 O, N 2 , CO, 


or H 2 ^. We are aware of only one study on acetone rad¬ 
ical fragmentation capable of also collecting neutral frag¬ 
ments. In specific, Ohkubo et ah— used a cryogenic STJ 
detector in a novel mass spectrometry approach for their 
investigations of 3 keV acetone radicals colliding with Xe 
atoms. They explained their data by two additional ETD 
channels producing solely neutral fragments: 

(CH3)2C0++A^CH3C0 + CH3+A+, (11) 

(CH3)2CO+-bH^CH3-bCH3-bCO-f A+. (12) 

However, their relative resolution of only AE/E r; 0.14 
was not sufficient to resolve molecular fragments differing 
only by the number of attached hydrogen atoms. 

Our results clearly show neutral fragments which do 
not originate from any of channels (nni) - (HU). The most 
interesting aspect of our data indicates that there is sig¬ 
nificant production of single C and O atoms as well as of 
the non-hydrogenated carbon molecules C 2 and C 3 . In 
the measurements of Ohkubo et al., due to their lower rel¬ 
ative energy resolution, any potential contribution from 
these fragments would be hidden in the peaks assigned 
as CH 3 , CO, or CH 3 CO. A full comparison between their 
work and ours certainly needs to take into account the 
differences between the two studies including collision en¬ 
ergy, target gas, and probably also internal excitation 
of the projectile ions. Nevertheless, the various neutral 
fragmentation channels are more clearly resolved in the 
present results. This demonstrates the wealth of new in¬ 
formation accessible by calorimetric detection of molec¬ 
ular fragments with high relative energy resolution. 

V. SUMMARY AND OUTLOOK 

We have studied the energy resolution function of an 
MMC detector for kinetic energy measurements of atomic 
and molecular projectiles in the energy range from ~ 13 
to 150 keV. For atomic projectiles we have demonstrated 
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FIG. 10. Same as Figure[Sl but for acetone radical ions (CH 3 ) 2 CO''’ at Sib = 150.45 keV (black line) and for neutral products 
from colliding these ions with residual gas (red curve). In both cases the Al-coated MMC absorber was used. The blue vertical 
lines mark the expected energies of the various neutral fragmentation products (see the text for a detailed explanation). The 
green arrow indicates the nominal ion beam energy Sib. 


that relative resolutions down to AE/E ^ 10“^ can be 
achieved. For C and H, the most relevant atomic pro¬ 
jectiles in organic chemistry, relative resolutions below 
10“^ were reached at all projectile energies studied. Ad¬ 
ditionally, we have studied the response of the MMC 
to molecular projectiles and demonstrated the capability 
of calorimetric mass spectrometry to identify molecular 
fragmentation products. 

We have also created a model of the micro-calorimeter 
energy resolution for molecular projectiles. The model 
was used to fit experimental data from and a 

very good qualitative and quantitative reproduction of 
the measured spectra was found. The quality of the 
model allowed us to resolve the contributions of the vari¬ 
ous energy-loss processes to the detector energy-response 
function. For molecules, the dominant distortion of the 
spectra originates from energy loss due to backscattering. 
Based on this finding we have experimentally demon¬ 
strated the strong influence of the absorber material on 
the detector energy-response function. Specifically, us¬ 
ing a thin aluminum coating on the absorber resulted 
in a strong reduction of the distortion from backscatter¬ 
ing. The validity of our model has been demonstrated for 
an MMC detector. Nevertheless, similar energy-loss pro¬ 
cesses are expected to occur in all types of cryo-detectors, 
including other micro-calorimeter types and STJs. Thus 
our findings on the MMC energy-response function are 
likely to be relevant for cryo-detectors in general. 

To the best of our knowledge, our work is the first 
experimental study of the energy-response function of 
cryo-calorimeter detectors for atomic and molecular pro¬ 
jectiles over the energy range covered. This model, which 
we have developed to estimate the energy-response func¬ 
tion, can be used by us and by others to extrapolate our 
findings to lower projectile energies and more complex 


molecular projectiles. In future studies we will investi¬ 
gate other coating materials in order to further improve 
the detector energy-response function. 

Implementation of micro-calorimeter detectors for 
fragmentation studies, such as the DR measurements at 
CSR proposed in Section HI requires not only a high en¬ 
ergy resolution, but also a large sensitive area, coinci¬ 
dent fragment detection capability, and position resolu¬ 
tion. To meet these requirements we have taken two 
different segmented MMC-design approaches. Segmen¬ 
tation provides the ability to detect coincident particles. 
The first approach consists of using large-area detector 
segments operated independently. Using a special con¬ 
figuration of thermal links between the absorber, sen¬ 
sor, and heat bath allows one to determine the incident 
position within each segment. This configuration has 
been realized and successfully tested on a circular pro¬ 
totype detector of 34 mm in diameter with 16 pie-like 
segments22i£i. The second approach involves our devel¬ 
opment of a large 4 kilo-pixel MMC detector with an 
active area of 44.8 x 44.8 mm^. The 64 x 64 particle ab¬ 
sorbing pixels making up this detector are coupled to a 
square array of paramagnetic temperature sensors. The 
sensor pickup coils are wired to provide for each row and 
each column of this array a summed temperature signal 
from all of the sensors in the respective row or column. 
This approach reduces the number of readout channels 
from 0{N‘^) to 0{2N), where N is the dimension of the 
square detector. Compared to the first approach, this 
pixelized detector will have a more uniform spatial reso¬ 
lution and higher multi-hit capability^^. 

Based on the results of our present MMC energy- 
response study and of the large-area position-sensitive 
MMC developments, briefly described above, we con¬ 
clude that MMC detectors are now a realistic option 
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for applications in atomic and molecular collision ex¬ 
periments and in mass spectrometry at keV energies. 
Our priority is to implement such a large detector in 
the recently commissioned electrostatic Cryogenic Stor¬ 
age Ring (CSR) at the Max Planck Institute for Nuclear 
Physics in Heidelberg, Germany^i^. The system will be 
used for electron-ion merged-beams studies as explained 
in Section m The < 10 K cryogenic environment of CSR 
facilitates the technical implementation of the large-area 
micro-calorimeter. The MMC energy resolution demon¬ 
strated here will be sufficient to efficiently detect and 
assign neutral products for DR of molecular ions stored 
in CSR. We expect achieve mass resolutions of 1 u for 
neutral DR products from parent molecular ions of mass 
^ 100 u and perhaps even heavier. The planned de¬ 
tector segmentation approach will allow for differentia¬ 
tion of the various DR fragments and their corresponding 
masses, thereby enabling the assignment of the various 
DR fragmentation channels (see Equation More¬ 

over, the position sensitivity connected with the detector 
segmentation will give us a measure on the kinetic energy 
released in the DR reaction for the various fragmentation 
channels. As a result, we will also be able to study the 
internal excitation of the product atoms and molecules 
can be investigated^. 

The micro-calorimeter detector, in combination with 
the CSR facility, will be a powerful tool for detailed 
studies of electron-capture reactions leading to neutral 
products (such as DR and ECD) in unprecedented de¬ 
tail, particularly for complex molecules. Moreover, we 
propose to extend molecular ion studies using micro¬ 
calorimeters to other reactions yielding neutral products, 
such as photodissociation, dissociative photodetachment, 
and collision-induced dissociation. While experimental 
techniques for studying some of these reactions exist 
(e.g.^S^^), they are usually unable to reliably distinguish 
and mass-assign multiple neutral fragments. The use of a 
segmented kinetic-energy-sensitive micro-calorimeter can 
overcome these limitations also for multi-channel frag¬ 
mentation reactions involving complex molecular ions, 
including species of biochemical interest. 
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